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SUMMAHf 


A method used by Tslen to derive similarity rules for hypersonic 
flows is utilized to derive von Kaiman’s similarity rules for transonic 
flows. A slight generalization is introduced by the inclusion of y, the 
ratio of specific heats, as a parameter. At the lower limit of the 
transonic region of flow the theory yields a formula for the critical 
stream Mach number's of a given family of symmetrical profiles. It Is 
further shown that this formula can also be obtained by means of the 
Prandt 1-Glaue rt small— pe rturbat i on method. Investigation of the behavior 
of the similarity parameter in the region where the thickness coefficient 
approaches zero and the critical stream Mach number approaches unity 
shows that it possesses a limiting value characteristic of the prescribed 
family of shapes. 


INTRODUCTION 


The rigorous solution of the subsonic flow of a compressible fluid 
past a prescribed closed body thus far has proved to be of insurmountable 
difficulty. As a consequence of this difficulty the emphasis has been 
placed on the establishment of a correspondence between the flow past a 
given body in an incompressible fluid and the same body in a compressible 
fluid. Among the best known results of this mode of attack are the 
Prandtl— Glauert rule and the von Karman-Ts i en velocity or pressure correc- 
tion factor — both based on some form of linearization of the f undame ntal 
nonlinear flow equations. None of the methods based on the linearization 
of the flow equations, however, can yield correct results in the transonic 
range where the flow is partly subsonic and partly supersonic. For this 
case a certain amount of the feature of nonlinearity of the flow equations 
must be retained in order to obtain useful and nontrivial results. In 
the present paper a detailed derivation is given of the transonic similarity 
rules recently given by von KarmAn (reference 1). 


FUNDAMENTAL EQUATIONS 


In plane steady flow the equation governing the flow of a nonviscous 
compressible fluid can be written in the form 
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In the derivation of this equation the pressure Is assumed to he a 
function of the density only. If further, the motion is irrotational, 
then 


= 0 ( 2 ) 


Here u and v are, respectively, the component velocities along the 
Cartesian x— and y-axes and c is the local velocity of sound given by 


c 2 = c 0 2 — J:(u 2 + v 2 ) (3) 

vhere cq is the velocity of sound at a stagnation point u = 0, v = 0 
and 7 is the ratio of specific heats at constant pressure and constant 
volume. 

Equation (l) is far too complicated to afford an insight into the 
properties of potential flow in the neighborhood of Mach number unity. 
The discussion is, therefore, restricted to the flow past a thin profile. 
Thus, at first v is assumed to be small in comparison with the sound 
velocity Cq. Equation (l) is then simplified to 


1 - 


y + 1 u 2 
• 1 u 2 


du 

5 ? 


dv 

5y 


c o 


(*) 


By the introduction of the sound velocity for which the local fluid 
and sound velocities are equal 


c* 


7 + 1 


'0 


and of the maximum possible fluid velocity 



(5) 


(6) 
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equation ( 4 ) "becomes 


1 - 


u 

C * 2 au + 


1 - 


U 


dx 


9max c 


dl = 0 

dy 


( 7 ) 


It is desirable to simplify this equation still further but yet to retain 
those features which yield nontrivial and useful results. Thus, 
if it is assumed that u is of the order of the critical velocity c* 

and if only teims up to the order 1 — J*- are retained, then equation (7) 

0* 


can be written as 


- (7 + l)(u - c* )~ + 0* = 0 

ox dy 


(8) 


Thus far the irrotationality condition, equation ( 2 ), has not been 
used. If, now, the undisturbed stream past a slender body is of velocity 
U slightly different from the velocity of sound and in the direction of 
the positive x— axis, then, according to equation (2) and the assumptions 
leading to equation (8), a velocity potential $ can be introduced with 


$ = c*x + (l — M*)cp 


( 9 ) 


where M* = — 
c* 


and (l — M*)cp is the disturbance-velocity potential. 


Then 


u = -& = c* + (1 - M*)i& 
dx dx 

v = M. » (1 _ M*)d$ 
dy dy 

and equation (8) becomes 


( 10 ) 


- (7 + 1)(1 - M*)d£ d§2> + c* d 2 ^ = 0 (11) 

dx dx 2 dy 2 

Equation ( 11) is a nonlinear simplified fom of the fundamental differential 
equation (l) and has been treated recently by von Kaman in connection with 
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similarity miles in two-dimensional transonic flow (reference 1). 

Equation (ll) * when expressed in hodograph variables., is of the type 
treated hy the Italian math e matician F. Tricon! some years ago (reference 2) 


DERIVATION OF SDCLAEITT CONDITIONS 


Recently* Tsien (reference 3) derived similarity rules for hypersonic 
flows where the fluid velocity is much larger than the velocity of sound. 

In the present paper the same procedure is employed to derive similarity 
rules for transonic flows where the fluid velocity is very nearly that of 
sound. 

According to the assumptions leading to equation (11)* it is implied 
that the solid body is thin, possesses no stagnation points, and the 
velocity of the fluid is everywhere in the neighborhood of the local 
velocity of sound. Now* suppose the profile of the obstacle to be 
symmetrical with respect to both the x— and y— axes and to possess cusps 
at both the leading and the trailing edges. Such profiles with uniform 
flow in the direction of the long axi s— of— symmetry x fulfill the 
assumptions leading to equation (ll). The flows past these profiles are 
said to be similar if the equation of motion (11) and the boundary 
conditions can be expressed in nondimensional variables in such a way 
that only a single constant factor is involved. Thus* if 2a is the chord 
and 2b is the maximum thickness of the body* then the following 
nondimensional variables are introduced: 


x = a| 

y = a(y + l) m t n ri 


> 


( 12 ) 


where t = — and m and n are exponents yet to be determined. It is 


clear that the nondimensional quantity involved is the thickness coeffi- 
cient t since this quantity determines the magnitude of the disturbance 
velocities. The exponents m and n are to be determined in such a way 
that the same constant factor appears in both the equation of motion 
the boundary conditions. 


The appropriate nondimensional form for the velocity potential <p is 


cp = ac*f(|, q) ( 13 ) 

By substitution from equations (12) and (13), the equation of motion (ll) 
becomes 
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(7 + l) aa+1 (l - M*)^ M & + & = 0 

di ^ dt f 


(1*0 


The boundary conditions at infinity require that the flow velocity 
he U. Hence, from the first of equations ( 10 ), 


g-— 


(at 00) 


The boundary condition at the surface of the slender body is 


(1 - M*)^$ = c*g^|)t (at y = 0 , -a < x ^ a) 

where g|— j describes the distribution of slope along the surface of the 
body. 

By means of equations ( 12 ) and (13), the boundary conditions can be 
written as : 


df 

as 

df 

dr) 


= -1 


r ( St 00 ) 


(7 + 1)^(1 - M *)^ 1 |£ = g( 6 ) (at q 


= 0, -1 ^ | < l) 


> (15) 


A comparison of the differential equation (14) and the boundary conditions, 
equations ( 15 ), shows that a single parameter is involved if 


2n - — (n + 1) or n - — i 

3 


2m + 1 = — m or m = — — 


that is. 
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(7 + 1) 1 /3( 1 _ M*)t~ 2 /3 = 


2K 


The undisturbed— stream Mach number ~ can be introduced in 

the following ways 

The general relation between M* and is 


Xl/2 


M* = 


1 + i M 2 
2 


(16) 


or 


Me 


= I 1 - (1 - M*)[2 - (1 - M*)' 


1 + 21—2(1 - M*) [2 - (1 - M*)] 


1/2 


Then, if powers of (l — M*) higher than the first are neglected. 


1 - M* = — ~( 1 - Moo) + . . . 


7 + 1 


Therefore 


(1 -Mj[(7 + l)t]“^ 3 = K 


(17) 


The results obtained thus far are such that by means of the substitution 
equations 


x = a| 

y = a [(7 + l)t] -1/ ^ 3 Ti 
<p = ac*f( I, tj) 

K = (1 - Ho) [(7 + l)t]“ 2/3 J 


(18) 


the differential equation for f(|, t|) and the boundary conditions became: 
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and 


_ rnr Sf 

' dt S|2 




(at oo ) 


*\ 




2K^ * g(!) (at 


rj = 0 , -1 £ | ^ 1 ) 


(19) 


( 20 ) 


The meaning of the similarity rule implied in the definition of the 
parameter K is the following: 

If a series of "bodies having the same distribution function g(| ) 
for the slope "but different thickness ratios t are placed in flowB of 
different undisturbed— stream Mach numbers and different va lue s of 7, 

such that the parameter E = 1 Z ft* , . . remains constant, then the 

[(7 + l)t] 2 ' 3 

flow patterns are similar in the sense that the same function f(|, tj) 
describes the flows. 


RESULTS DERIVED FROM TEE SIMILARITY RULE 
Pressure Coefficient 


In the case of a uniform flow past a fired boundary, the pressure 
coefficient is defined as 


C P>H» = 






where and are, respectively, the pressure and density in the 

undisturbed stream and the static pressure p in the fluid is given by 


P = P 
00 



2 00 


( 


u^ + v^ 
U 2 


JL. 

7-1 


( 21 ) 
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Then 





r 


*\ 

2 

—1 4- 

l 7 ~ hi 2 ( u2 + - y2 i) 

7-1 L 

11 2" 
7M» 


2 ' 00 \ ) 







( 22 ) 


By means of equations (10) and (l8) , if powers of (1 — M^,) higher than 
the first are neglected, the following result is obtained 


p „ k i-n.( 

P ’"“ - - — —TT( 


i + H) 
/ 


or 


2/3 


:p ' M "“i^(V !i) i A P(s ’ n;K) 


(23) 


where P(|, r\; K) depends on the form of the solution f (| , q) for 
the particular family of profiles treated. 


Lift Coefficient 


The lift l of the body 1 b given by 


l 



dx 


By a similar procedure, as in the derivation of equation (23), the lift 
coefficient is given by 


c l 


l m JL l(k) 

|pJJ 2 (2a) Me, 2 (777)173 


(2k) 


where 


Pi 

L(K) = P(|, 0; K) d| 

J“1 
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and where for an extremely thin straight— line profile the thickness 
coefficient t has been replaced by the angle of attack a. 


The pressure drag d 


d 

Hence the drag coefficient 


Drag Coefficient 

of the body is given by the following expression-, 
0 a 

= 2 J-a ** 

is given by 


c d - 


d 

^U 2 (2a) 



£/3 

(7 + D 1/3 


E(E) 


(25) 


where 


■u 

D(K) - 2 g(|)P($, 0; K) d| 

0-1 


ADDITIONAL COKIDERATIOie 


The results derived in the present paper apply to two-dimensional 
near-sonic flows past thin shapes. Such flows have been calculated 
for a family of symmetrical shapes with cusped leading and trailing 
edges (reference 4 ) and for a family of elliptic cylinders (reference 5). 
These calculations are valid at least up to the critical stream Mach 
number M^. The critical Mach number may be considered from two points 
of view. First,, it may be considered to denote the lower limit of 
a mixed subsonic— supersonic flow, that is, where the imbedded supersonic 
region is simply the point of maximum fluid velocity at the surface of 
the solid. From this point of view, according to equation (17), a 
critical value of the similarity parameter K can he defined. Thus 


Kc 


r 


1 - Her 

( 7 + l) t] 


( 26 ) 


This equation can also be written in a form that yields the critical stream 
Mach number for a given family of shapes ; that is. 
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Mcr = 1-Kc r [(r + l)t] 2/3 (27) 

Second* the critical— stream Mach number may he considered to denote the 
upper li mi t of the purely subsonic range of speeds. This point of view 
suggests a derivation of equation (27) by means of the Prandtl-Glauert 
small— perturbation method (reference 6 ). The procedure is as follows; 

The relation between the local and the undisturbed— stream Mach 
number* within the approximation of the small— disturbance theory* is 
given by 


M 2 = 


1 + 


2 ir ( 1 + 


M. e ) 


( 28 ) 


where u 1 is the disturbance velocity and U is the undisturbed— stream 
velocity. 

By definition M = 1 for Mj,, = Mq^. and, since C-p = —2— 
approximately * equation (28) becomes 


1 =Mcr 2 



( 29 ) 


The relation between Cp^j/[ cr and the pressure coefficient Cp,o 
of the incompressible fluid is given by the Prandtl-Glauert rule 


C p,M cr 


C ^° 

Vl - Mcr 2 


Hence equation (29) becomes 


, 2x3/ 2 

C 1 - Mcr9 c 


( 30 ) 


As seen from equation (29)* 1 — M,^ 2 is of the first order in the sma l 1 

perturbation 3 ^- and* accordingly* equation (30) includes terms of higher 

order. This fact can be seen by rewriting equation (30) in the following 
foim: 
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(1 - M or 2 ) 3/2 

' L ~ [l “ (l - Mcr 2 ) ] [l " ffj (1 ~ **cr 2 )_ 




Then, to within the lowest order in the small perturbation, 

l-^r 2 = (- 2 T J «P.o ) 2/3 

or, approximately, 

«or = 1 -f(- 2 4 Jc P' 0 ) 2/3 (31) 

It is quite easy to show the connection between equation (31) and 
equation (27) . Thus, for an incompressible fluid, the pressure coefficient 
is given by 


C P,0 


x-s§ 

u 2 


(32) 


where q is the magnitude of the fluid velocity at any point in the 
field of flew. 

In the case of the family of shapes of reference k s the maximum 
velocity at the surface is given by 


^ = 1 4 + « . # 

U 2 

Hence 


'P*0 


- -3t + 


and equation (31) becomes 



( 33 ) 
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In the case of the family of elliptic cylinders of reference 5, 
the maximum velocity at the surface is given ty 

§ = i + t 


Hence 


'P,0 


= -2t + 


and equation (31) becomes 


“cr = l-£[(7 + l)t] 2/3 (3*) 

An examination of equations (33) and ( 34 ) suggests that 
possesses a limiting value; that is. 




lim 1 ~ Her 
t ■* 0 [ (y + l)t] 2/3 

Mcr-*1 1 J 


(35) 


It is noteworthy that the value of Spends on the family of 

shapes , although in the limit t -* 0 the profile in every case is a 
straight-line segnent. The numerical values of (K cr ) 11wi shown in 

equations {33) and (34) represent, however, only the effect of the 
Prandt 1 -Glaue rt of first-order term in the power-series development in t 
of the max fauna velocity at the solid surface® It is rather surprising 
that the higher-order terms involving the higher powers of t also 
contribute to the value of ( K cr) n Tn ‘’ 1116 P roce4ur ® is singly to replace 

the maximum velocity by the critical speed; that is. 
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where a^, ag, a^, . . . depend on the given family of shapes and 
involve only 7 and 


Then as suggested by equation 
2Kcr [(? + l)t J equation (36) 


(26), when 1 — M^ 2 is replaced by 
is identically satisfied for t 0, 


with each additional higher-order teim contributing to the value of K or » 

In this fashion there are obtained, successively, linear, quadratic, 
cubic, and higher-order equations for the determination of |^K cr | . |im . 


The foregoing considerations have been applied to the family of 
symmetrical shapes of reference 1 and to the family of elliptic cylinders 
of reference 5. The values of for values of t different from 

zero were obtained from equation (26) with the required values of 
determined by means of the theoretical results of references 1 and 5* 

The limiting values of were obtained from equations corresponding 

to equation (36). Tables I and II show the results of these calculations 
for air and for Freon-12. Note the approach to the limiting values of 

as t approaches zero and approaches unity. The first column 

of values of Kc r shows the effect of the first-order or Prandt 1-Glaue rt 

term; the second and third columns show, respectively, the effects of the 
second-order and third-order terms in the thickness coefficient t„ The 
successive values of Kqp for the various values of the thickness 
coefficient t indicate good convergence. 


Langley Memorial Aeronautical laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va„, December 9* 19^7 
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TABLE I 

VALDES OF K cr FOE THE FAMILY OF SYMMETRICAL SHAPES OF HEFEEENCE 4 
[For air, 7 = 1,4; for Freon-12, 7 = I.I36] 


t 

M cr 

K„„ 

cr 

Approximation 

Approximati on 

First 

Second 

Third 

First 

Second 

Third 

Air 

0. 

,002 

.005 

.010 

.020 

.040 

.060 

.080 

,100 

1.000 

.982 

.967 

.947 

.918 

.872 

.837 

.807 

.781 

1.000 

.980 

.963 

.942 

.909 

.859 

.818 

.784 

.753 

1.000 

• 979 
.962 
.940 
.906 
.854 
.812 
.774 
.743 

0.655 

.647 

.638 

.632 

.624 

.609 

.595 

,581 

.568 

0.717 

.710 

.705 

.700 

.689 

.675 

.664 

.650 

.639 

0.738 

.732 

.726 

.719 

.709 

.696 

.685 

.678 

.665 

Freon-12 

0 

.002 

.005 

.010 

.020 

.040 

.060 

,080 

.100 

1.000 

.983 

.969 

.951 

.924 

.881 

,848 

.819 

.794 

1.000 
'• .981 
.966 
.946 
.915 
.867 
.830 

.797 

.769 

1.000 

.980 

.964 

.944 

.912 

.863 

.823 

.787 

.756 

0.655 
.647 
.642 
.636 
.627 
.612 
.600 
.588 
• 576 

0.717 
.713 
.709 
.704 
.694 
. 682 
.670 
.658 
.647 

O.738 
.735 
.730 
.725 
.716 
,706 
.698 
.691 
.68 2 
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TABLE II 

VALUES OF K cr FOE THE FAMIIX OF ELLIPTIC CYLINDERS OF REFERENCE 5 
air, 7 = 1 . 4 ; for Frean- 12 , 7 = 1,136] 



^or 

® cr 

Approximati on 

Approximati on 

First 

Second 

Third 

First Second 

Third 


Air 


0 

1.000 

1,000 

1,000 

0.500 

0.537 

0.556 

,002 

.986 

.985 

.984 

.494 

.530 

.550 

,005 

.974 

*972 

.971 

,490 

.526 

.544 

,010 

,960 

.957 

*955 

.487 

.522 

.537 

,020 

.936 

*932 

*930 

.481 

*515 

.526 

.040 

.901 

.894 

.892 

,470 

*505 

.515 

,060 

. 873 

,863 

.861 

,462 

,497 

.507 

,080 

.848 

*837 

.833 

.455 

.490 

.500 

.100 

.827 

.813 

,809 

.447 

.484 

.494 

Freon-12 

0 

1,000 

1.000 

1,000 

0,500 

0.537 

0.556 

,002 

.987 

• 986 

.985 

.497 

.534 

.553 

,005 

*976 

.974 

*973 

*495 

.532 

.550 

,010 

,962 

.959 

.958 

,491 

,528 

.543 

.020 

.941 

. 936 

*935 

,484 

.520 

.534 

,040 

.908 

.900 

.899 

.474 

.509 

.522 

,060 

.881 

,872 

.871 

.469 

,504 

.513 

,080 

.858 

.847 

.845 

.462 

.496 

.506 

.100 

.838 

,826 

.822 

*455 

.488 

.498 
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